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ABSTRACT 

The spectral and magnetic properties and variability of the B2Vnp emission-line magneto- 
sphere star HR7355 were analyzed. The object rotates at almost 90% of the critical value, 
meaning it is a magnetic star for which oblateness and gravity darkening effects cannot be 
ignored any longer. A detailed modeling of the photospheric parameters indicate that the star 
is significantly cooler than suggested by the B2 spectral type, with R e ^ = 17 500 K atypi- 
cally cool for a star with a helium enriched surface. The spectroscopic variability of helium 
and metal lines due to the photospheric abundance pattern is far more complex than a largely 
dipolar, oblique magnetic field of about 11 to 12 kG may suggest. Doppler imaging shows 
that globally the most He enriched areas coincide with the magnetic poles and metal enriched 
areas with the magnetic equator. While most of the stellar surface is helium enriched with 
respect to the solar value, some isolated patches are depleted. The stellar wind in the circum- 
stellar environment is governed by the magnetic field, i.e. the stellar magnetosphere is rigidly 
corotating with the star. The magnetosphere of HR 7355 is similar to the well known a Ori E: 
the gas trapped in the magneto spheric clouds is fairly dense, and at the limit to being optically 
thick in the hydrogen emission. Apart from a different magnetic obliquity, HR7355 and the 
more recently identified HR5907 have virtually identical stellar and magnetic parameters. 

Key words: Stars: individual: HR7355 — Stars: early-type, magnetic field, chemically pe- 
culiar 



1 INTRODUCTION 

The helium- strong star HR7355 (Rivinius et al. 2008 ) has recently 
been found to host a magnetic field with several kilogauss longi- 
tudinal strength (< B z >) (Rivi nius et al.||2010[ |Oksala et al] 
12010] ). The rotational period of P = 0.52 d puts it among the 
shortest period non-degenerate magnetic stars known, while its 
v sin i « 310 km s _1 is the highest for any known non-degenerate 
magnetic star. 

For a long time, rapid rotation (the term "rapid rotation" used 
here as being at a fraction of well above 50% of break-up velocity) 
and the presence of a strong magnetic field have been regarded as 
mutually exclusive for early type stars, and so the sheer existence of 
this object poses a number of puzzling questions. Most important 
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is the stellar age vs. the expected spin down time. Age estimates for 
HR7355 range from 15 to 25Myr (Mi kulasek et al.|2010| ), much 
longer than the same authors give for the spin-down timescale, and 
hence clearly problematic for fossil field explanations. However, 
the gravity of logg = 3.95 found by Oksala et al. (2010) would 
be more supportive of a younger star. In any case, HR7355 is a 
star for which effects like gravity darkening and oblate deforma- 
tion cannot be ignored. This means traditional analysis methods to 
derive stellar parameters will, at best, give uncertain results, and at 
worst misleading ones. 

Intriguing is the presence of photospheric abundance nonuni- 
formities, supporting the result that rotational mixing even in rapid 
rotators should be inhibited at such field strengths (Zahn 201 1) . 
Not only are there abundance variations across the stellar surface 
in HR7355, but the amplitude of the Hel equivalent width (EW) 
variations is much larger than in the similar, though less rapidly 
rotating, Bp star a Ori E. 

For the following work we will briefly review the rotational 
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ephemeris in Sect. then introduce the new observations in 
Sect. [2] The first aim of this work is to provide stellar parameters 
taking into account the rapid stellar rotation, i.e. rotational defor- 
mation and gravity darkening, in Sect. [3] The spectroscopic vari- 
ations are then described and analyzed in detail in Sect. [4] as are 
the measured properties of the magnetic field in Sect. [5] and of the 
magnetosphere in Sect. [6] The results are discussed in Sect. [7] and 
finally the conclusions are drawn in Sect. [8] 



Table 1. Summary of UVES observations 

Run date range # of spectra 

(MJD) civil (UT) blue/red arm 

A 54923 - 5495 1 2009-04-02 to 04-30 29/33 
B 55026-55062 2009-07-14 to 08-18 55/59 
C 55077 - 55083 2009-09-03 to 09-09 20/20 



1.1 Ephemeris 



There are several ephemerides and period determinations, out of 
which only the one by Riviniu s et al.|p008] > is clearly obsolete. In 
order to converge on one, we combined the data used by Rivinius 
(|2010| ) with the photometric data published by Mikulasek 



et al. 



et al 



et al 



(2010), and the photometric data kindly provided by Oksala 



pOK)] ). While the period by |Rivinius eTal] pOK)] ) did not 



sort well the combined photometric data when phased with this 
period, the 2008 photometry by |Oksala et~aT] and the 2008/2009 
spectroscopy (Rivinius et al. and this work) do have their minima 
at acommon phase. As a result, we adopt here the epoch given by 
Rivinius et al] ( |2010| ) T (MJD) = 54 940.33, which is the mid- 
date of an occultation of the star by the magneto spheric material, 
seen in the B aimer line spectroscopy, and as such very well defined. 
However, the better period proved to be the value published by Ok- 
sala et al.[ P — 0.5214404(6) d, who could rely on a time baseline 
twice as long. The error of the period is such that for the purposes 
of this work it can be fully neglected, and we treat the period as a 
fixed value. 



2 OBSERVATIONS 

Apart from archival and literature data, which will be referenced 
where used, this work is based on high-resolution echelle spec- 
tra obtained in 2009 with UVES ( |Dekker et al.|2000) at the 8.2m 
Kueyen telescope on Cerro Paranal. The instrument was used in the 
dichroic 2 437/760 setting, which gives a blue spectrum from about 
375 to 498 nm and a red spectrum from about 570 to 950 nm, with 
a small gap at 760 nm. The slit width was 0.8", giving a resolving 
power of about R — 50 000 over the entire spectrum. The UVES 
observations were obtained in three service mode observing runs, 
in April, July, and September 2009 (see Table[T]for the date ranges). 

Exposure times were between 120 seconds in April and 30 
seconds in July to September, the latter ones being repeated four 
times in a row. Since the period is short we did not average shorter 
exposures, but treated them as taken at distinct phases. A few of the 
blue spectra taken with t e x P = 120 sec were overexposed in better- 
than-average seeing conditions, thus we have more red than blue 
spectra. In total, we have 104 blue and 112 red spectra. The typi- 
cal signal to noise ratio (S/N) of a single spectrum is about 275 in 
the blue region and 290 in the red. However, for unknown reasons, 
high exposure levels in UVES, even if still in the fully linear regime 
of the detector, tend to cause some order-merging and global nor- 
malization problems, i.e. high S/N data obtained with UVES often 
have more residual echelle order "wiggles" than data with a lower 
exposure level. This prompts for local re-normalization of the spec- 
tra for analysis of any particular line. In particular, broad spectral 
features, like Balmer-line wings, are to be treated with care. Our 
analysis of those features is cross-checked on the low resolution 
FORS1 long-slit data of |Rivinius etahlpOTO] ). 
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Figure 1. The averaged UVES line profile for C11A 4267 and the best fitting 
model for this particular equivalent width, with T e ff = 18 500 K, and the 
final best model, taking into account additional conmstraints as discussed 
in the text, with T eff = 17 500 K, both with vsini = 310 km s -1 and 
shifted to v sys = 7kms _1 . 



3 STELLAR PARAMETERS 

In order to determine the stellar parameters, we made use of spec- 
tral synthesis, both for line strength and profiles, as well as for ab- 
solute fluxes. The code we used is the third, re-programmed version 
of |Townsendl s {1997) BRUCE and KYLIE suite, which in the fol- 
lowing we will call "Bruce3". 

Apart from computational issues, and the implementation of 
perturbations like pulsation, which is not used here, the main im- 
provement is the treatment of limb-darkening: the code uses either 
NLTE atmosphere models by |Lanz & Hubeny (2007]), or |Kuruczj s 
(T992) LTE (ATLAS 9) atmospheres, from which the synthetic 
spectra were computed by SYNSPEC ( [Hubeny & Lanz|201 1| >, with 
wavelength-dependent limb-darkening coefficients, i.e. including 
the effects of limb-darkening on the spectral lines, ranging from 88 
to 750 nm. Several characteristics of the line profiles, most notably 
strengths of Cll and Sill lines in the blue range, were better repro- 
duced with the LTE atmospheres rather than with the NLTE ones 
(but still using an NLTE spectral profile code). The spectrum of 
HR7355 is, therefore, modeled based on an ATLAS9-SYNSPEC 
grid. 

For line-profile comparison, the synthetic spectra, although in 
principle already "perfectly" normalized (i.e. according to the com- 
puted continuum), were re-normalized with the same procedure as 
the observed spectra, i.e. with a series of smoothly joining spline 
segments, and using the same wavelengths to sample the contin- 
uum. 
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3.1 Projected rotational velocity and radius 

In a first step, the averaged profile of CllA 4267 was used to obtain 
the projected rotational velocity by comparing the profile to the grid 
of Bruce3 models. For any reasonable set of parameters, v sin i = 
310±5 km s _1 approximated the width of the observed profile best. 
We note that this approach completely neglects macroturbulence 
or other potential broadening mechanisms. However, concerning 
macroturbulence this is a reasonable approach for a main sequence 
star, and significant magnetic broadening is not expected for the 
measured field. The value is in good agreement with Oksala et al. 
pOTO] ) (see Fig. IT). 

The systemic velocity was determined during this measure- 
ment as well, to v sys — 7 db 2 km s _1 . We found no indication for a 
variation of this value, hence we can exclude a companion close or 
massive enough to significantly affect the magnetosphere, unless it 
would have an orbit in the plane of the sky, i.e. significantly mis- 
aligned with the rotation, seen equatorially (see below). Assuming 
rigid rotation of the circumstellar environment, the derived v sin i 
in combination with the rotational period of P = 0.5214404 d the 
equatorial radius then becomes i?*, eq sini = 3.19 Rq. Assuming 
that a rapidly rotating star is constrained by five independent pa- 
rameters, usually v eq , i, T e ff , M*, R*, eq , this leaves only three of 
them to be determined, namely the mass, the effective temperature, 
and the inclination. 



3.2 Effective temperature 

The interpretation of the effective temperature T e ff in a rapidly ro- 
tating star, i.e. in a star where gravity darkening and oblateness have 
become significant, is not straightforward. So we note that Bruce3 
uses T e ff in the sense that it is the uniform temperature, that a black- 
body of the same surface area would need to have, to have the same 
luminosity as the actual gravity darkened star, integrated over the 
full solid angle, while one has to keep in mind that observation- 
ally derived T e ff are usually solid-angle dependent quantities and 
only assumed to be valid over the entire sphere. The value for T eq 
and Tpoi in Table [2] are the temperatures, as given by the gravity 
darkening relation (Colli ns|1963| >, at the stellar equator and pole, 
respectively. The plane-parallel input atmospheres are interpolated 
to these values. 

Again the line of CllA4267 is most suited for constraining 
T e ff, at least when relying on NLTE line modeling, because its 
strength strongly depends on T e ff, while the modeled values of 
i or M* turn out to hardly have an effect on either strength nor 
shape. We note that systematic uncertainties of the behaviour of 
CllA 4267 may affect the value, however, within the frame of using 
the SYNSPEC input model spectra, as computed above, the value 
is well constrained. For this step, solar abundance values were as- 
sumed. The best fit with a Bruce3 compute SED is obtained for 
T e ff = 18 500 K on the UVES mean spectrum (see Fig. [TJ, and 
T e ff = 18 250 K on the FORS1 mean spectrum. The fits for other 
spectral lines are reasonable as well, but since those do vary more 
strongly with i and M* than CllA4267 they have less power in 
constraining T e ff. 

An alternative, independent way to constrain T e ff is to use the 
International Ultraviolet Explorer (IUE, Bog gess et al.|1 978) flux- 
calibrated spectra and visual photometry, de-redden them with an 
appropriate extinction law, and then fit the model fluxes. This not 
only gives the temperature by fitting the slope of the UV spectrum, 
but also the distance by fitting the absolute flux as measured above 
the Earth's atmosphere vs. the absolute flux as modeled at a unit 



distance from the star. The latter will be used in a later step, for 
now we perform the exercise only on the slope to obtain T e ff . Two 
spectra exist, both taken in high resolution with the short wave- 
length arm (SWP camera). As for both the large aperture was used, 
they can be considered as spectrophotometric, and are, in fact, in 
good agreement with each other. Both spectra were obtained from 
the IUE archive, co-added and finally median filtered to reduce the 
noise. For the photometry we used the Stromgren values given by 
|Gr0nbech & Qlsen| ( ^97oT ) and the flux calibration of |Gray| ( [T998] ). 

The value for the color excess E^ B _ V ) — 0.08 mag was taken 
from Papaj & Krelowski (1992). However, they adopt a standard B2 
star as template, having an intrinsic color of (B — V)o = —0.21, 
which due to the He- strong nature of the object is not appropriate. 
There is consensus that HR7355 is in fact cooler than a typical B2 
star, n amely about 17 to 18 kK (Mik ulasek et al.|2010||Oksala et al.| 
2010} , and thus an intrinsic color of (B — V)o = —0.18 might be 
a better estimate, so that the excess needs to be adjusted to a lower 
value of E^ B _ V ) — 0.05 mag. The effect of the magnetic field on 
the colors at the temperatures, metallicities, and field strenghts in 
question is entirely neglibible ( Kochukhov et al.|2005| >. 

To actually de -redden the IUE data, we use the extinction 
law in the parametric form given by |Cardelli et al. ( 1989|, as- 
suming a standard Galactic value of Rv = E^ B _ V )/Av — 3.1. 
An independent determination of Rv would require spectropho- 
tometry of the 2200 A extinction bump, but unfortunately no IUE 
long- wavelength data were taken. We note, however, that changing 
Rv does not modify the conclusions significantly, and in partic- 
ular does not provide any help in the distance problem presented 
below. We find that the slope of the UV data between 1250 and 
2000 A is fitted best with T eff - 17 000 K. Since this seems a bit 
low compared to the line determinations, we tried to fit a model 
for T e ff = 18 000 K as well, but found this is only possible if we 
adopt a higher reddening of E^ B _ V ) — 0.08 mag, which is argued 
above to be implausible. In addition, it should be kept in mind that 
the star is chemically peculiar, after all, so that an analysis based 
on the assumption of solar abundance might show some systematic 
effects. 

Since the measurements differ of T e ff somewhat, at this point 
we go ahead with both options T e ff = 17 000 K and T e ff = 
18 000 K, and will decide at the end of the procedure which to adopt 
for the final model. 

Although the cool temperature derived here seems incompati- 
ble with the spectral classification of B2, usually assumed to be at 
about 22 000 K, it was mentioned above that the He- strong nature of 
the star is responsible for the discrepancy. Early B stars are mainly 
classified on account of their Hel lines, and since the helium lines 
are much stronger than they would be for a star with solar abun- 
dances, classifications of B-type He-strong stars are systematically 
biased towards the spectral type with the strongest Hel lines, which 
is B2. 

3.3 Inclination and stellar mass 

In terms of Bruce3 modeling, since all other relevant quantities, 
like period and v sin z, are well constrained, the inclination effec- 
tively sets the model's equatorial radius (since it converts vsini 
to the rotation velocity) and the geometric projection onto the line 
of sight. As the equatorial radius is thus constrained, the model- 
ing parameter range is chosen to include masses between 6 and 
8Mq, considering stellar evolution: lower mass stars at the main 
sequence are not large enough for the required equatorial radius, 
neither are higher mass stars small enough. The considered range 
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Figure 2. The Bruce3 model (see Table[2} vs. the observed data: Upper panel: SED fit for UV IUE short-wavelength spectrophotometry and Stromgren visual 
photometry (large red dots) vs. modeled fluxes. Middle panel: Line profile fit for IUE spectrum (red) vs. modeled UV line profiles (black). Lower panel: Line 
profile fit for averaged UVES (red) and FORS1 (blue) spectra vs. modeled visual line profiles (black). Note that in this step solar abundances are still used, 
hence the mismatch for the helium profiles. This discrepancy is solved later by the Doppler imaging. 



is sufficiently broad to to rely on details of the evoolution, it is, af- 
ter all, just the range over which a grid is computed. The possible 
inferred range of effective temperature supports this mass range. 

In terms of modeling, the stellar mass only scales the polar ra- 
dius and the amount of gravity darkening. However, both quantities 
have only subtle effects on the modeled spectral lines, and results 
are not really conclusive. However, in one respect the model does 
vary strongly with these two parameters, namely in the total lumi- 
nosity and in the flux seen by the observer, since they scale the size 
of the stellar surface through the radius and the surface tempera- 
ture distribution. Therefore, the traditional approach, determining 
the distance with the modeled flux, can be inverted: The revised 
Hipparcos parallax is sufficiently precise to do so. The Hipparcos 
distance, 273 ± 26 pc, is used to define a range for model fluxes; 
namely "near" (247 pc), "mean" (273 pc), and "far" (299 pc) dis- 
tances. 

We find the models for % — 80° do not have sufficient flux; 
they remain well under the faintest flux limit, even for the "near" 
distance. For T e ff = 18 000 K, models for i = 70° are just in 
agreement with the "near" distance flux, and only models with i — 
60° are compatible with the flux curve for the "mean" distance. For 
T e ff = 17 000 K, even models with i = 60° are just compatible 
with the flux curve for the "near" distance. The "far" flux is too 



luminous to be reproduced by any of the computed parameter sets. 
The ultimately- adopted inclination is, therefore, i = 60° for which 
we assume a ±10° uncertainty. 

Although this still does not constrain the mass in a straight- 
forward observational way, all the other parameters are sufficiently 
well-known to leave only a narrow range of physically acceptable 
masses. This is because we can, at least to first order, expect that 
the luminosity of a star of a given mass does not change with its 
rotational velocity. 

The luminosity computed by Bruce3 modeling, however, is 
not constrained by the physical processes inside the stellar interior; 
it is just given by the temperature distribution and the stellar sur- 
face area, both of which can be chosen independently in Bruce3. 
The stellar evolution requirement that mass and luminosity are re- 
lated is an independent constraint, taken into account only in a later 
step (and as well for setting a reasonable grid range above). This 
means the Bruce3 -computed luminosity is almost invariant with 
model mass. Instead, it increases with decreasing inclination, in 
such a way that all T e ff = 18 000 K, i = 60° models have a lumi- 
nosity of about L* = 1120 L , T eff = 17 000 K, % = 60° models 
have a luminosity of about L* = 1000 Lq. Such a luminosity is 
only compatible with a mass of about 6 ± 0.5 Mq. Unfortunately, 
both models give a very similar distance to match the absolute flux, 
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Figure 4. Phased variations of photospheric helium and metal lines, represented as dynamic spectra. Top panels: original profiles, the dotted lines mark the 
cuts applied for black and white in the grey-scale panels. Bottom panels: phased dynamical spectra. These panels were constructed with 32 or 64 phase bins 
(64: SimA4553 and SnA4816, HeiA6678, HeiA4713) depending on the complexity of the variations prompting for a trade-off of higher sampling vs. the 
noise level. For a detailed description of the variations see Sect.|4] 



around 240 pc, which is only somewhat closer than the Hipparcos 
distance of 273 ± 26 pc. 

So the only remaining question, left open above, is whether 
to adopt 17 000 or 18 000 K. While the model for 18 000K agrees 
better with the spectroscopic results, and has a more plausible lu- 
minosity, it would require a somewhat higher reddening than sug- 
gested by E( B - v) ■ On the other hand, the model for 17 000 K pro- 
vides a better fit for the spectral energy distribution, yet its lumi- 
nosity would be lower by 10 to 20 % than expected even for a star 
of 6 Mq only. However, given the relatively low extinction this is 
a rather robust determination. For the final model, therefore, we 
adopt 17 500 K, but unfortunately have to accept a rather generous 
error margin of 1000 K. 

The final parameters are summarized in Table [2] the absolute 
flux fit for these parameters is shown in Fig. [2] The figure illustrates 
that these parameters give good fits to both absolute flux and line 
profiles, except of course the variable helium lines, as well as some 



of the UV lines. Concerning the uncertainty intervals given in Ta- 
ble^ we note that only the one given in the uppermost part (input 
parameters) were constrained traditionally, i.e. by judging a grid of 
model vs. observational data. The intervals given for the derived 
parameters, however, were obtained by propagating these input pa- 
rameter uncertainties through the model code, and thus they do not 
necessarily satisfy the usual statistical relations (in particular, they 
are non-orthogonal to each other). Furthermore, since the quanti- 
ties have been computed, we give more significant digits than the 
interval would warrant, to keep the values of the derived parameters 
consistent with the input parameters. 



4 PHOTOSPHERIC ABUNDANCE PATTERN 

Figure [5] shows previously published data and equivalent widths 
measured in the UVES data, phased with the ephemeris given in 
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Table 2. Stellar parameters. Only errors given in the "input" part indicate 
goodness of fit as discussed in the text, errors in the "derived" part are 
strictly the intervals derived from propagating the errors of the input pa- 
rameters. Since the uncertainty of the period is insignificant wrt. that of the 
other quantities we treat it as fixed. 



Bruce3 input parameters 



Stellar mass 
Inclination 
Effective temperature 
v sinz 

^Prot 







6±0.5M 
60 ± 10° 
17 500 ± 1000 K 
310 ± Skins" 1 
0.5214404 d 



Bruce3 derived parameters 

358 kms- 1 ± 10% 
0.89 ± 7% 



veq 

W = Veq/V cr it 

Polar radius 
Polar temperature 
Polar gravity (log) 
Equatorial radius 
Equatorial temperature 
Equatorial gravity (log) 
Luminosity 



3.06 R ±5% 
19 751 K ±7% 
4.25 ± 11% 
3.69 R ± 10% 
15 740 K ±9% 
3.85 ± 13% 
972 L ±28% 



Adopted spectrophotometric parameters 
Excess E^ B _y^ 0.065 mag 

Distance 236 pc 



Sect. The uppermost panel shows the available photometric 
data. A constant offset on magnitude has been added to the ad- 
ditional data sets, taken in different photometric systems, for the 
purpose of display together with the Hipparcos data. We note that 
photometric data from different bands may not have the same am- 
plitude, as well as shifts in phase (expressed here a linear phase 0) 
vs. each other, but the accuracy of the available data is not sufficient 
to determine if this is the case here. 

In any case, the data by |Oksala efal]p010| ) and the UVES 
spectra were taken only one year apart, and the period is known 
well enough to be sure now about the phasing of the spectroscopic 
variations vs. the photometric: When Ha is in occupation, there is 
a photometric minimum. However, since the photometric minima 
are broad and the phases of occultations and Hel EW minima differ 
only by A<j> « 0.07, the following can also be said: When the Hel 
absorption is weakest, there is a photometric minimum. Whether 
the photometric variations are due to the occultations, as suggested 
by[Tbwnsend (2008), or due to the photospheric abundance pattern, 
as suggested by Mikulase k et al.|(2010) , will require careful mod- 
eling, and possibly further multicolor photometric observations to 
distinguish geometric from intrinsic effects. 



4.1 Line strength variations 

The Ha equivalent widths are fully dominated by variations of 
the circumstellar environment; there is no indication of an intrinsic 
photospheric variation of this or any other hydrogen line strength 
(see Fig. [7]). The vertical dotted lines in the four upper panels of 
Fig.[3]mark the central points of the occultations, i.e. when the cir- 
cumstellar lobes are centrally in front of the stellar disk (at v sys ). 
This is at = 0, actually fixed by the definition of the ephemeris, 
and at 4> = 0.46. 

The third and fourth panels in Fig. [3] show the equivalent 
width variations of the HelAA4713 and 4388 lines. The double 
wave character of the Hel equivalent width curve is very clear. Al- 
though the strongest EW points in both half-cycles are of a similar 
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Figure 3. Phased observational data. Top panel: Photometric magnitudes 
from Hipparcos (green), Oksala et al. ( 2010 binned to 64 phase bins, blue 
open circles), and AS AS data (red triangles, binned to 16 phase bins) from 
Mikulasek et al. (2010). Second, third, and fourth panels: Equivalent width 
measurements of Ha, HeiAA4388, and 4713 in mA, in the UVES data 
(black circles), the 2008 FORS1 spectra jRivinius et al.|2010| purp le trian- 
gles), and the 1999 and 2004 FEROS spectra (Rivinius et al.|2008] red and 
blue squares, respectively). Lowermost panel: Magnetic data from FORS1 
( Rivinius et al.|2010] see there a version of this plot with less data), and 
ESPaDOns magnetic data (open circles) in gauss and a sinusoidal fit to the 
FORS1 measurements. The vertical lines denote special phases discussed 
in the text. 



level, the weakest EW values differ by a factor of two. The dashed 
lines mark the estimated times of maximal He-line strength at about 
= 0.25 and 0.68. 



4.2 Line profile variations 

The most obvious photospheric variation is that of the Hel lines. 
The variations are very strong, relative to the line strength much 
stronger than in the prototype a Ori E: e.g. the EW of Hel A 4713 
varies between 250 and 550 mA in HR7355 vs. between 420 and 
570 mA in a Ori E (measured in our own and archival data for both 
stars). Figure [3] shows a factor of two and four between the EWs of 
the weak and strong states of Hel lines. Judging the EW variations 
in the figure at face value, they form a simple double wave pattern 
that would be well compatible with just two enhanced polar regions 
and a depleted belt. A detailed look at the spectra, however, shows 
that the variations are not that easily explained. 

The forbidden transition of Hel at 4045 A (Fig. [4] upper left) 
shows the simplest pattern; presumably because it is the weakest 
line of Hel only the strongest features are seen: Two He-enhanced 
regions are visible, which, according to the phasing with the mag- 
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netic curve in Fig. [3] are consistent with being located at about the 
magnetic poles, as is typical for He-strong stars (see e.g. |Groote| 
|& Hunger|[T982l |Reiners et al.|2000| for a Ori E). These two re- 
gions are also the dominant ones in other helium lines. They cross 
^sys = = 7kms _1 at 0o — 0.19 and 0o — 0.65. 

However, the stronger lines show additional features. In 
HelA 4713 and 6678 (Fig. [4] middle two panels in upper row) three 
more, albeit weak, features are seen: They cross zero velocity at 
0o = 0.02, 0.46, and 0.88. The position of the first two agree with 
the occultations (see Sect.[6|. However, the time these features take 
to propagate across the profile (formally an inverse acceleration, 
in dynamical spectra directly visible as "slope" of a propagating 
feature), show this is not the case: The time to cross the entire 
range of 2 v sin i for undoubted occultation features, in all lines, 
is much faster, as seen at 0o = 0.46 in the panels for Nil A 4631 
and SilllA 4553 in Fig. [4] This is so because the magnetosphere 
is empty (or at least far less dense) closer to the star than about 
2 stellar radii, as seen by the lack of Ha emission at velocities 
< 2 x v sin i at quadrature, and the crossing time of a corotating 
feature scales inversely with the distance from the star. 

In the fairly many Nel lines detected in the spectrum, like 
NelA6717 shown in Fig. [4] (upper right panel), there are two fea- 
tures crossing the profile, at 0o = 0.20 and 0o = 0.64. These 
are almost the same phases as for the strong helium features. Nel, 
therefore, is a metallic species that is probably enriched at the mag- 
netic poles, rather than at the magnetic equator, as one usually ob- 
serves for metals in He- strong stars. 

Unfortunately, looking at more metal lines does not make the 
picture any simpler: While in Nil A 4631 only two features are seen 
(0 O = 0.35 and 0.93) four each can be identified in CllA 4267 and 
Sim 4553 (<£ = 0.08, 0.32, 0.49, and 0.95), and possibly as many 
as six are seen in SllA 4816 (<£ = 0.02, 0.13, 0.32, 0.49, 0.60, and 
0.73). 



4.3 Doppler imaging 

We applied the Doppler imaging (DI) technique to interpret phase- 
dependent behavior of some of the variable spectral lines in terms 
of an inhomogeneous distribution of chemical elements. Due to 
the strong rotational Doppler broadening of its spectra, HR 7355 
is a very challenging target for DI. In fact, it has the largest v sin i 
among the early-type stars to which abundance DI has been ap- 
plied. Due to excessive blending and shallow line profiles, only 
CllA4267 and strong He lines, such as HelA4713, are suitable for 
modeling. We used these two spectral features to derive respective 
abundance distributions with the updated version of the magnetic 
DI code described by Piskunov & Kochukhov (2002). 

Since HR7355 has a very strong magnetic field, it has to be 
taken into account in abundance mapping. Available Stokes I ob- 
servations do not allow the derivation of the magnetic field geom- 
etry simultaneously with chemical spots. We therefore adopted the 
11.6 kG dipolar magnetic field model derived in the next section. 
On the other hand, gravity darkening and oblateness of the star 
were not included in the modeling. As a consequence, we might 
expect spurious pole-to-equator abundance gradients for the lines 
especially sensitive to effective temperature. For the prupose of DI, 
stellar parameters, i.e. inclination and v sin i, were adopted as de- 
rived above form spectroscopic modeling. 

The results of the DI analysis of CllA 4267 are illustrated in 
Fig. [5] The surface distribution of carbon required to reproduce 
variability of this line features significant meridional gradient, pos- 
sibly due to unaccounted gravity darkening. In addition, there is 



about 1 dex abundance contrast at high latitudes, dominated by two 
spots seen at rotational phases 0.25 and 0.80. 

Helium is overabundant in the atmosphere of HR7355 and 
shows a high-contrast surface distribution resulting in dramatic 
variability of Hel lines. Since concentration of He is comparable 
to that of H over a significant part of the stellar surface, not only 
is the observed line profile variability a result of the abundance de- 
pendence of the He line absorption coefficient, but is is also a con- 
sequence of continuum instensity variations and, possibly, lateral 
variations of the atmospheric structure. All these effects have to be 
accounted for in order to obtain an accurate absolute scale and mag- 
nitude of the He abundance variations. We carried out He mapping 
with the help of a new magnetic DI code, Invers13 ( Kochukhov 
|et al.|2012) , which uses a grid of model atmospheres to treat self- 
consistently the impact of abundance variations on the line profiles. 
Our present calculations are based on the LLmodels ( Shulyak et al. 
12004} grid computed for T eff = 17500 K, \ogg = 4.0 and the 
helium abundance from log (He/ H) = —3.0 to log(He/H) = 1.9. 

The results of He abundance inversion are presented in Fig. [6] 
The overall contrast of the He map is nearly 4 dex. The main He 
overabundance features are located at the equator and below. To 
reproduce a narrow bump traveling across the He line profile at 
phases 0.95-0.05, the code requires a latitudinally extended spot 
where He is strongly underabundant relative to the Sun. 



5 THE MAGNETIC FIELD 

While the high v sin i tells that the inclination cannot be small, 
some further constraints on both the inclination i and the oblique- 
ness of the magnetic field vs. the rotational axis f3 can be obtained 
from the observed properties of the longitudinal magnetic field. 

The magnetic field curve is not centered around zero, i.e. 
B m [ n ^ -5 max , as noted in the original papers ( |Oksala et al.| 
|2010|[Rivinius et al.|2010| ). If a dipole is assumed, this means nei- 
ther i nor p can exactly be 90°, although they might be close: If 
the inclination was 90°, the rotational projection of the two mag- 
netic poles must be the same, regardless of /3. The same is true for 
P — 90°, since then the magnetic poles would be on the rotational 
equator. 

Although the field, for reasons discussed in Sect. [7] is most 
likely not exactly dipolar, the variation is sufficiently sinusoidal to 
Preston] fl971| > to the magnetic variation 



apply the equations by 



( |2010| >. This yields r — B m j n /B max — 



given by |Rivinius et al. 

cos(/3 + i)/cos(P - i) = -0.78. If the values of |Oksala et~aT. 
(2010) are used, r = —0.80 is obtained, but this changes the val- 
ues of P and i (see below) only by a few degrees. This value of r 
means that P + i must be smaller than or equal to 141°. If the sum 
is equal to 141°, then i = P = 70.5°. If i ^ P, we can say that 
on the one hand, i cannot be much smaller than 70° because of the 
large v sin z, on the other hand i cannot be much larger either, since 
then P would have to be very small, which would make it very hard 
to explain the strong mean longitudinal field, unless the dipole field 
strength reaches several dozen kilogauss or more at the poles. We 
conclude that both i and P are most likely between 60 and 80°, in 
such a way that ft + i is about 130 to 140°. This limit to the incli- 
nation is in nice agreement with the spectroscopic determination, 

which is 60 ± 10°. 

U sing |Sch warzschild s (1950) in the form given by |Leone] 
|et al.| ( |2000| ), and assuming a pure dipole with i — 60° and 
P = 75°, the field strength at the magnetic poles would then be 
about 11.6kG. 
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The observational indications for a non-dipole field and a non- 
sinusoidal field curve are subtle, so that the above values, summa- 
rized in Table [3] should at least be a good approximation. 



6 THE CIRCUMSTELLAR MAGNETOSPHERE 

The circumstellar magnetosphere is filled with magnetically 
trapped gas (see e.g. Figs. 1 and 8 of Town send et al.|[2007] for 
an illustration of the concept and definitions). The trapped gas, ac- 
cording to Townsend et al. occupies a warped plane between the 
magnetic and rotational equator, effectively the points along each 
field line where the centrifugal acceleration is highest, as long as 
this is larger than the local gravity. However, this plane is not uni- 
formly filled; the intersections of the magnetic and rotational equa- 
torial planes are preferred, forming two geometrically relatively 
thin lobes ( Townsend 2008 ). This gas is seen in emission in hydro- 
gen lines when in quadrature and in absorption in hydrogen and a 
few other lines when passing through the line of sight. Species and 
lines in which the occultations are present include Hi, OlA 8446 
and the 777x triplet, Sill A 6347, Nil A 4631 (see respective panel 



Table 3. Magnetic parameters. The inclination 7 
Table[2] 



: 60 ± 10° is taken from 



Field parameters 
(3 75 ± 10° 

Polar B dipole 11.6kG 

Magnetospheric lobe parameters 
Inner radius 2 R* 

Outer radius 4 R* 

Particle density io 12 - 3 ±°- 5 C m~ 3 
Plane tilt < 22°/ cos % 

Lobe separation 165 ± 5° 



in Fig. [4] at 4> — 0.46 and v = Okms 1 ), and many Fell lines 
(FellA 4949 is visible in the SilllA 4553 panel in Fig.|4] at = 0.46 
and v = -200 km s" 1 ). 
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Figure 6. Same as Fig.[5]for helium distribution reconstructed using the HeiA 4713 line. The He abundance is given in logarithmic units relative to the Sun. 



6.1 Geometry 

The two occultations occur at phases <j) — 0.0 and <j) — 0.46 
(Fig. [7}. The ph enomena associated with the passing of the mag- 
netospheric lobes in front of the star provide information on the 
circumstellar geometry. The lobe absorption is a relatively nar- 
row feature, of a FWHM of only about 250 km s -1 in Pai4. This 
value is similar in other lines. This corresponds to about 40 % of 
2 x (vsini). Since the circumstellar gas is locked to the magne- 
tosphere, it is in co-rotation with the stellar surface. For this rea- 
son any projected radial velocity interval is unambiguously mapped 
onto a strip parallel to the projected rotational axis, and the width 
of such an interval translates directly to geometry: A 90° tilt would 
produce zero width (as only one velocity would be obscured at 
any given time), while a 0° tilt would obscure all velocities across 
2 x (v sin i) simultaneously. Hence, the projection of the tilt of the 
plane of the lobes against the rotational axis must then be around 
the arcus tangens of the simultaneous velocity dispersion of the cir- 
cumstellar absorption in units of 2 x (vsini), which in this case 
is arctan 0.4 = 22°. This actually is an upper limit, if a non-zero 
scale height for the lobes was assumed, the inclination had to be 
even lower. Note, however, that this tilt is not just a projection of 
90° — P, since the plane of the lobes is somewhat warped. Never- 
theless, it is probably close to this value. 



An equatorial feature in bound corotation, traversing the line 
of sight at photospheric level, would take half a cycle to do so 
(A(p = 0.5). With increasing distance, the crossing time, i.e. the 
fraction of the orbit within a line of sight impact parameter < 1, 
decreases quadratic ally. The passing of the lobe through the line of 
sight towards the star is much faster than Acj) = 0.5: from — v sin i 
to +v sin i takes the feature about Acj) = 0.13. i.e. about one fourth 
of the crossing time as in case material would be present down to 
the stellar surface. ThereforeCombining this and the quadratic de- 
crease with distance, it can be concluded that there is no absorbing 
material directly above the star, but sufficient density to absorb is 
found only at and above a radius of 2 R* . 

Outside the occultation phases, the emitting material is seen 
next to the stellar disk, and again since it is in magnetically bound 
corotation there is a linear relation between rotational velocity and 
distance from the stellar surface, meaning material at twice the ra- 
dius has twice the rotational velocity. Thus, when the circumstellar 
material is in quadrature, the observed velocities give additional 
geometry information: The emission has an inner edge at about 
600kms _1 , which again points to an empty region inside 2R*. 
The outer edge of the emission is at a velocity of about 4 x v sin i 
in Ha and 3 x v sin i in Pai4, which gives outer limits for the lobe 
emission of about 4 R* and 3 R*, respectively. 

There is also a "skew" in the appearance of the emitting lobes, 
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Figure 7. Observational variations of the circumstellar magnetosphere. The upper row shows the actual profiles overplotted, in the lower row the phased 
variations are shown. For Hen and H/3 the residuals after subtracting the Bruce3 model are used for the phased plot; the residual level in the line center is zero 
to better than 0.01. Emission is clearly seen in the residuals, in fact up to HS. Paschen lines show emission as well; shown here is Pai4 (note the different 
scales of the abscissae). 128 phase bins were used for all panels in order to sample the fast variations around the occultation phases in sufficient detail. 



Table 4. Equivalent width of the circumstellar emission in the B aimer lines 
and the computed B aimer decrements for two phases and the red and blue 
emission peaks separately. 







Equivalent widths [A] 


B aimer deer. 






Ha 


H/3 


H 7 


D34 D54 


Blue, 4 


) = 0.17 


-1.37 


-0.42 


-0.17 


1.16 0.59 


Red, 4 


> = 0.17 


-0.94 


-0.25 


-0.13 


1.33 0.76 


Blue, 4 


) = 0.78 


-1.00 


-0.33 


-0.14 


1.07 0.62 


Red, 4 


) = 0.78 


-1.17 


-0.26 


-0.08 


1.59 0.45 



well seen in Fig. ^] Tracing the maximum emission levels belong- 
ing to each lobe, it is found that each lobe varies in a sinusoidal 
shape, i.e. the cardinal points of conjunctions and quadratures are 
separated by a quarter of a cycle, but the two lobes are separated 
by less than half a cycle, i.e. geometrically by less than 180°. That 
difference is the same as the one that causes the second occultation 
to be at 4> = 0.46 instead of at = 0.5. 

6.2 Density 

The theoretical line profiles are a very good approximation of the 
mean photospheric profiles, as seen in Fig.|2]for H7 and Fig.[7]for 
H/3 and Ha. Therefore, the residuals faithfully represent the pure 
circumstellar emission and it is possible to apply nebular diagnos- 
tics and calculate B aimer decrements. The B aimer decrements D34 
and D54 were measured separately for the blue- and redshifted 
emission peaks at phases 0.17 and 0.78, when the emitting mate- 
rial is seen next to the star (see Table j4jl. For computi ng the decre- 
ments, the same method as used by Stefl et al. (2003) for the clas- 
sical Be star 28 CMa was adopted. The stellar continuum fluxes to 
correct for the normalization were those modeled by Bruce3, so 
that /*(Ha)//*(H/?) = 0.35 and /*(H 7 )//*(H/3) = 1.46. The 
measured values were compared to the theoretical computations by 



Williams & Shipman (1988), which were derived for an isother- 
mal pure hydrogen accretion disk of 10 000 K, that is optically thin 
in the continuum. The B aimer decrements depend only weakly on 
temperature in such a hot plasma, so that even if the circumstellar 
material of HR7355 is hotter than 10000K the derived densities 
still provide useful guidelines for modeling. While the values for 
D54 would be in agreement with logarithmic particle densities be- 
tween 11.7 and 12.5 per cm 3 , the logarithmic particle densities in- 
ferred from D34 are somewhat higher, at 12.2 to 12.8 per cm 3 . In 
any case, these values are close to the optically thick limit, above 
which the decrements become independent of density. 



7 DISCUSSION 

7.1 HR 7355 as a rapid rotator 

HR7355 was the first magnetic star discovered to rotate so rapidly 
that rotational effects like oblateness and gravity darkening should 
be taken into account. A second example, HR5907, has been pub- 
lished by (Gru nhut et al.|2012) . Due to the relative proximity of 
HR7355 not only was it possible to obtain high quality data, but 
the public archives contain data for the star as well. 

Using all these data, the five parameters required to describe a 
rapidly rotating star could be derived (Table [2}. Of these five, two 
stand out: Firstly, the effective temperature of just 17 500 K. This 
is fairly low for a spectral classification of B2. We note, however, 
that B2 is the spectral subtype at which Hel lines are strongest, and 
so the classification for a star with helium overabundance will be 
biased towards that spectral subtype with respect to its temperature. 
Although the Doppler imaging indicates that not the entire surface 
is overabundant in helium; with patches of both over- and under- 
abundance existing, the stellar surface on average is overabundant 
in helium, justifying its classification as He- strong in earlier papers 
( |Oksala et al.125T5||Rivi5us et al.|2010] >. 
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The second unusual parameter is the rotation rate, at almost 
90% of the critical value (even the previous record holder for 
rapid rotation among non-degenerate magnetic OB stars, a Ori E 
barely reached 50%). This means that no significant rotational spin- 
down can have occurred for HR7355 yet. Unless the field is self- 
generated and has only existed since recent times, this leaves two 
possibilities: Either the star is so young that it had no time to brake 
yet, or the mass loss is so weak that it does not remove any signif- 
icant amount of angular momentum. While the latter would have 
been unlikely for a star of 22 000 K, i.e. corresponding to a normal 
B2 star, it is much more plausible for a star of only 17 500 K. For 
the very similar HR5907, Gru nhut et al . | \20 1 2] ) obtain a spin-down 
timescale of 8Myr. This is still less than the 15 to 25Myr sug- 
gested by Mikulasek et al. (2010). However, for an age of less than 
8 Myr our stellar parameters, i.e. gravity, mass, polar radius, and 
luminosity, agree with evolutionary models. Even at a distance of 
only 236 pc the conclusion by Mikulasek et al. (2010), namely that 
HR7355 is probably not a member of Scorpius-Centaurus, stands. 

7.2 Magnetic field geometry 

While for a photospheric parameter analysis the tools exist to take 
rotation into account, such tools for the magnetic analysis are only 
being developed and have not yet reached the parameter space re- 
quired here, simply because the inclusion of such strong effects 
was not necessary to date. The most rapid rotation taken into ac- 
count for magnetic analysis so far was by Khalack (2005), limited 
to Rpoie/Reci < 0.95. In the following discussion, the effects of 
the rapid rotation are speculated about, but we stress firm conclu- 
sions will require more sophisticated theoretical work to be done. 

A sinusoidal variation of the longitudinal magnetic field is the 
hallmark of a centered dipole field tilted with respect to the ro- 
tation axis. However, the high-quality spectroscopic data suggest 
that the field topology may be more complex than this. With the 
circumstellar clouds situated at the intersections between magnetic 
and rotational equators, the cloud transits should be coincident with 
B z — 0. The measured magnetic curve, when fitted with a sine, has 
an offset of +0.3 kG. The magnetic nulls are expected to occur at 
= and 4> = 0.54 (dotted lines in lowermost panel of Fig. [3j. 
However, the cloud transits are observed to occur at = 0.0 and 
<j) = 0.46 (dotted lines in upper panels of Fig. [3}. The maximal 
strength of the Hel lines normally coincides with the magnetic 
poles pointing towards the observer. This would be expected for 
(j) — 0.28 and 0.78 in case of a sinusoidal variation (dashed lines in 
lowermost panel of Fig.[3j, but spectroscopically, i.e. estimated by 
maximal He-enhancement, occurs at = 0.25 and 0.68 (dashed 
lines in upper panels of Fig. [3}. 

Offsets as discussed in the previous poaragraph could theo- 
retically be explained in case of a de-centered magnetic dipole. 
The two misaligned cloud crossings would require an offset of the 
dipole along the magnetic polar axis. Such an offset could also ac- 
count for the constant term of 0.3 kG in the field curve, although 
this alone could also be explained by the projection topology of 
an aligned field. The observed timing of the magnetic poles fac- 
ing towards the observer requires an offset perpendicular to this 
axis. These offsets would show in a non- sinusoidal variability curve 
of the magnetic field. In any case, however, off-center dipoles, or 
any combination of low order multipole components, can only be a 
simplified approximation to the actual magnetic field, as magnetic 
Doppler imaging has shown that more complex and small scale 
structures are usually present even for stars with largerly sinusoidal 
variations in the longitudinal field ( Kochu khov et al.|2011| >. 



The line profile variations are surprisingly complex, with at 
least four, possibly as many as six subfeatures identifiable in some 
lines. It is surprising that a relatively strong and ordered field as 
in HR 7355 gives rise to so many clearly distinct zones of helium 
and metal abundances. One possibility to be explored could be an 
interplay of a more simple type of abundance variation pattern with 
the co-latitude dependent line formation in a rapidly rotating star 
due to the gravity darkening. 



7.3 The magnetosphere 

The circumstellar magnetosphere is filled with dense gas, suffi- 
ciently dense to make the B aimer emission optically thick. Since 
the gas needs to be ionized to be trapped by the magnetic field, 
one might wonder whether in the central regions of this dense gas 
this is still the case. The occultation pattern observed in other line 
transitions can be used for guidance to answer this question: The 
signature of the occultation of the central, most dense part of the 
magnetosphere is seen in lines of Fell, Nil, Sill etc. It is not seen, 
however, in the lower excitation/ionization lines of Call or Nal. 
This means that the temperature of the magnetosphere is still of the 
order of 10 000 K, where hydrogen is still far more than sufficiently 
ionized to be locked to the magnetic field. 

7.4 Comparison to a Ori E and HR 5907 

The prototype of an emission line magnetosphere star is a Ori E. 
It is, in fact, fairly similar to HR7355: It is He-strong, with simi- 
lar polar field strength and i = 75° and f3 = 55°. The two main 
differences concerning parameters are the higher effective tempera- 
ture (23 000 K) and slower rotation (less than 50% critical) of a Ori 
E (see e.g. Oksala et al.|2012[ for the most recent parameter deter- 
mination). As a consequence of geometry and field strength, the cir- 
cumstellar variations are comparable, except for the lower v sin i. 
Fig. 5 of Oksal a et al.| is extremely similar to our Fig. [7] As well, 
occultation crossings can be seen in metal ions of a Ori E, as in 
Fig. 6 of that work, FelllA5127 at = 0, vs. our Fig.[4j NllA4631 
at (j> — 0.46. A morphological difference, for which it is not clear 
how it is associated with the parameters, is that the equivalent width 
changes due to the surface chemical abundance pattern are of much 
lower amplitude in a Ori E than in HR 7355. 

Shortly after HR7355 was confirmed to be a magnetic star 
( Oksala et al. 2010 ; Rivinius et al. 2010 ), a very similar object with 
an even shorter period was found, HR 5907 ( Grun hut et al.|201l] 
2012 ). Mass, effective temperature, and inclination of HR 7355 and 
HR 5907 are within each other's error bars (see Table[2|, and the pe- 
riod as well hardly differs with 0.521 d vs. 0.508 d. The equatorial 
velocity of HR 5907 is somewhat slower, and the equatorial radius 
smaller than for HR 7355, so that HR 5907 is less luminous, but the 
average photospheric spectra of HR7355 and HR5907 are virtu- 
ally identical. Also the dipolar field strength for both is of the order 
of 10 kG, and the gas density of the magnetosphere is comparable 
as well. 

The only strong difference is the obliquity of the field, being 
much smaller for HR5907 (almost aligned, /3 = 7°), which, how- 
ever, has striking effects on the variability. The amplitude of the 
line profile variability is much smaller in HR 5907, since the mag- 
netic equatorial belt is visible to the observer all the time, and the 
contributions from the pole only vary slightly, while for HR7355 
magnetic pole-on and equator-on views alternate. Nevertheless the 
actual profile variability is still much more complex than in slower 
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rotators, even if the equivalent width amplitude is very small. As 
pointed out by Townsend (2008), the low obliquity also leads to 
only one occultation per rotational cycle, while HR7355 has two. 
As well due to this the emission variability of HR 5907 is not a 
clear double-helix pattern as for HR7355. 



8 CONCLUSIONS 

We have shown that HR7355 is an early B type magnetosphere 
star. In principle similar to the well known a Ori E, it stands out 
as the non-degenerate magnetic star rotating most closely to the 
critical limit, at about 90%, followed by HR5907. The magnetic 
field of HR7355 has a polar strength of about 11 to 12 kG, with a 
morphology close to, but as suggested by the asymmetries in the 
equivalent width curves probably not exactly dipolar, and with a 
fairly large obliquity of 75°, seen at an inclination of about 60°. 

The stellar parameters point to a relatively young star with an 
effective temperature typical for a mid B rather than an early B 
type star, namely 17 500 K. As the typical mass-loss rate for such 
a temperature is significantly less than for a B2 type star, the spin- 
down timescale is several Myr, and hence the age and rapid rotation 
are not inconsistent with the presence of a fossil magnetic field. 

The stellar surface shows strong abundance variations, which 
have been Doppler imaged for helium and carbon in this work. Al- 
though some small surface areas in the Doppler imaging have a 
helium abundance lower than the solar value, helium enriched ar- 
eas dominate. Indeed, compared to a model with typical solar or 
B-star abundances, the Hel lines are stronger in the time averaged 
profile. This is also the explanation for the spectral classification as 
type B2 in spite of the effective temperature of only 17 500 K. 

As recently as ten years ago, a Ori E was the only B-type mag- 
netic star with an emission line magnetosphere, and 1 OriC was 
the only such O star. In the last few years this situation has changed, 
and HR7355 is not only a new member of the growing class of 
magnetosphere emission line stars, but has doubled the parameter 
space for these objects in the rotational dimension. Density esti- 
mates for the magneto spheric region are similar in all three B stars 
mentioned (<r Ori E, HR 5907, HR 7355). 
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